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Abstract: 

The complete matrix elements for e"'"e~ —>■ 4/ and e"'"e~ — 4/7 are calculated in the 
Electroweak Standard Model for polarized massless fermions. The matrix elements for 
all final states are reduced to a few compact generic functions. Monte Carlo generators 
for e^e~ — >■ 4/ and e"'"e~ — >■ 4/7 are constructed. We compare different treatments of 
the finite widths of the electroweak gauge bosons; in particular, we include a scheme with 
a complex gauge-boson mass that obeys all Ward identities. The detailed discussion of 
numerical results comprises integrated cross sections as well as photon-energy distributions 
for all different final states. 
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1 Introduction 

When exceeding tlie W-pair production threshold, the LEP collider started a new era 
in the verification of the Electroweak Standard Model (SM): the study of the properties 
of the W boson and of its interactions. While the most important process at LEP2 in 
this respect is certainly e^e~ W"^W~ — > 4/, many other reactions have now become 
accessible. Besides the 4-fermion-production processes, including single W-boson produc- 
tion, single Z-boson production, or Z-boson-pair production, LEP2 and especially a future 
linear collider allow us to investigate another class of processes, namely e"'"e~ —>■ 4/7. 

The physical interest in the processes e"'"e~ — * 4/7 is twofold. First of all, they are 
an important building block for the radiative corrections to e"''e~ —>■ 4/, and their effect 
must be taken into account in order to get precise predictions for the observables that 
are used for the measurement of the W-boson mass and the triple-gauge-boson couplings. 
On the other hand, those processes themselves involve interesting physics. They include, 
in particular, triple-gauge-boson-production processes such as W"*'W~7, ZZ7, or Z77 
production and can therefore be used to obtain information on the quartic gauge-boson 
couplings '-ffWW, '-fZWW, and j'jZZ. While only a few events of this kind are expected 
at LEP2, these studies can be performed in more detail at future linear e"'"e~ colliders [ 

i- 

Important contributions to e+e —>■ 4/ and e"'"e — > 4/7 arise from subprocesses with 
two resonant gauge bosons. These subprocesses do not only contain interesting physics, 
such as the quartic gauge-boson self-interactions, but also dominate the cross sections 
in those regions of phase space where the invariant masses of certain combinations of 
final-state particles are close to the mass of the corresponding nearly resonant gauge 
bosons. Therefore, it is sometimes a reasonable approach to consider only the resonant 
contributions. In the most naive approximation, the produced gauge bosons are treated 
as stable particles. Most of the existing calculations for triple-gauge-boson production 
have been performed in this way, such as many calculations for e"'"e~ —>■ W"''W~7 [ ^. 

In an improved approach, the so-called pole expansion [ ^, ^, the resonances are 
treated exactly, i.e. the decay of the produced gauge bosons is taken into account, but the 
matrix elements are expanded about the poles of the resonances. Once the non-resonant 
diagrams have been left out, the expansion is in fact mandatory in order to retain gauge 
invariance. If only the leading terms in the expansion are kept, this approach is known as 
the pole approximation, or double-pole approximation in the presence of two resonances. 
The accuracy of the (double-)pole approximation is, at best, of the order of Ty/My, 
where My and Fy are the mass and the width of the relevant gauge bosons, and thus 
typically at the level of several per cent. Consequently, this approach is only reasonable 
if the experimental accuracy is correspondingly low. The error estimate of Ty/My is too 
optimistic in situations in which scales of order Ty, or smaller, are involved. This is, in 
particular, the case in threshold regions or if photons with energies of the order of Ty 
are emitted from the resonant gauge bosons. On the other hand, the quality of the pole 
approximation can be improved by applying appropriate invariant-mass cuts. 

Note that the double-pole approximation is particularly well-suited for the calculation 
of the radiative corrections to gauge-boson pair production, since the error resulting from 
the double-pole approximation is suppressed by an additional factor a/vr in this case. In 
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this approximation the corrections can be classified into factorizable and non-factorizable 
corrections The virtual factorizable corrections can be composed of the known 

corrections to on-shell W-pair production [ Q and on-shell W decay [ 0, and the non- 
factorizable corrections have recently been calculated [ |^. In contrast to the virtual 
factorizable corrections, the real factorizable corrections cannot be simply taken over from 
the on-shell processes. In the case of real photon emission the definition of the double- 
pole approximation is non-trivial if the energy of the final-state photon is of the order 
of Fv, since the resonances before and after photon emission are not well separated in 
phase space. A possible double-pole approximation of the 0{a) corrections to four-lepton 
production has been discussed in Ref . [ |^ . 

In our calculation of e"'"e~ 4/7 we do not use the double-pole approximation since 
the full calculation is feasable with reasonable effort. Unlike the double-pole approxima- 
tion, the full calculation is valid for arbitrary processes in the set e+e" — >■ 4/7. Moreover, 
when using the exact results for the real corrections in an 0{a) calculation of 4 fermion- 
production processes, the reliability of possible approximations can be tested. 

Some results for e"'"e~ 4/7 with an observable photon already exist in the litera- 



ture. In Refs. [ |T0], ^ the contributions to the matrix elements involving two resonant 



W bosons have been calculated and implemented into a Monte Carlo generator. This 



generator has been extended to include collinear bremsstrahlung [ [T^ and used to discuss 
the effect of hard photons at LEP2 [ [l^. The complete cross section for the process 
e"'"e~ ude~UeJ has been discussed in Ref. [0]. In Ref. [0, the complete matrix 
elements for the processes e"'"e~ — > 4/7 have been calculated using an iterative numerical 
algorithm without referring to Feynman diagrams. We are, however, interested in explicit 
analytical results on the amplitudes for various reasons. In particular, we want to have 
full control over the implementation of the finite width of the virtual vector bosons and to 
select single diagrams, such as the doubly-resonant ones. So far no results for e"'"e~ — > 4/7 
with e"'"e~ pairs in the final state have been published. 

In order to perform the calculation as efficient as possible we have reduced all processes 
to a small number of generic contributions. For e"'"e~ —>■ 4/, the calculation is similar to 



the one in Ref. [|T6[, and the generic contributions correspond to individual Feynman 
diagrams. In the case of e"'"e~ — 4/7 we have combined groups of diagrams in such a 
way that the resulting generic contributions can be classified in the same way as those for 
e+e~ 4/. As a consequence, the generic contributions are individually gauge-invariant 
with respect to the external photon. The number and the complexity of diagrams in the 
generic contributions for e"'"e~ —>■ 4/7 has been reduced by using a non-linear gauge-fixing 
condition for the W-boson field [ . In this way, many cancellations between diagrams are 
avoided, without any further algebraic manipulations. Finally, for the helicity amplitudes 
corresponding to the generic contributions concise results have been obtained by using 



the Weyl-van der Waerden formalism (see Ref. [ |T8| and references therein). 

After the matrix elements have been calculated, the finite widths of the resonant 
particles have to be introduced. We have done this in different ways and compared the 
different treatments for e"'"e~ 4/ and e"''e~ 4/7. In particular, we have discussed a 
"complex-mass scheme" , which preserves all Ward identities and is still rather simple to 
apply. 
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The matrix elements to e^e~ 4/ and e"'"e~ 4/7 exhibit a complex peaking 
behaviour owing to propagators of massless particles and Breit-Wigner resonances, so 
that the integration over the 8- and 11-dimensional phase spaces, respectively, is not 
straightforward. In order to obtain numerically stable results, we adopt the multi-channel 



integration method [ |I6|, |T9| and reduce the Monte Carlo error by the adaptive weight 
optimization procedure described in Ref . [ . In the multi-channel approach, we define a 
suitable mapping of random numbers into phase space variables for each arising propaga- 
tor structure. These variables are generated according to distributions that approximate 
this specific peaking behavior of the integrand. For e^e~ — >■ 4/ and e"'"e~ — 4/7 we iden- 
tify up to 128 and 928 channels, respectively, which necessitates an efficient and generic 
procedure for the phase-space generation. We wrote two independent Monte Carlo pro- 
grams following the general strategy outlined above. They differ in the realization of a 
generic procedure for the construction of the phase-space generators. 

The paper is organized as follows: in Section ^ we describe the calculation of the 
helicity amplitudes for e"'"e~ —>■ 4/ and e"'"e~ —>■ 4/7 and list the complete results. The 
Monte Carlo programs are described in Section |^. In Section ^ the numerical results are 
discussed, and Section ^ contains a summary and an outlook. 

2 Analytical results 

2.1 Notation and conventions 

We consider reactions of the types 

e'^iP+, + e~ip-, (T-) fiih, (Ti) + /2(A;2, cxa) + fsih, (T3) + f^ih, ai), (2.1) 
e'^iP+, + e"(p_, o-_) ^ fi{ki, ai) + f2{k2, ^2) + fsih, 0-3) + U{h, a^) + 7(^5, A). (2.2) 

The arguments label the momenta p±, ki and helicities = ±1/2, A = ±1 of the 
corresponding particles. We often use only the signs to denote the helicities. The fermion 
masses are neglected everywhere. 



For the Feynman rules we use the conventions of Ref. [ ^ . In particular, all fields 
and momenta are incoming. It is convenient to use a non-linear gauge-fixing term [ ^ of 
the form 

2 



fix 



+ ie(A^ - — Z'^)iy+ - iM, 



--(9^Z,-Mzx)^--(9M,)^ (2.3) 

where 0^ and x are the would-be Goldstone bosons of the and Z fields, respectively. 
With this choice, the ip^W^A vertices vanish, and the bosonic couplings that are relevant 
for e"'"e~ — >■ 4/7 read 



V 1 




—\egvww [di^pik- — k^)^ — 2g^ykv,p + '^g^pky^u] 
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= -2iegvww g/iuOpa, (2.4) 

with V = A, Z, and the couphng factors 

gAww = 1; 9zww = — —■ (2.5) 

Note that the gauge-boson propagators have the same simple form as in the 't Hooft- 
Feynman gauge, i.e. they are proportional to the metric tensor Qf^iy. This gauge choice 
eliminates some diagrams and simplifies others owing to the simpler structure of the 
photon-gauge-boson couplings. 

The vector-boson-fermion-fermion couplings have the usual form 




ie7ME^?y/,/,^-' (2.6) 



where uj± = (1 ±75)72. The corresponding coupling factors read 



s P 1 



where Qi and j = ±1/2 denote the relative charge and the weak isospin of the fermion 
/i, respectively, and // is the weak-isospin partner of The colour factor of a fermion 
/, is denoted by N^, i.e. A/'fepton = 1 and N^^^^^ = 3. 

2.2 Classification of final states for e+e~ — 4/ 

The final states for e^e~ 4/ have already been classified in Refs. [ |2^. We 



introduce a classification that is very close to the one of Refs. [ . It is based on 

the production mechanism, i.e. whether the reactions proceed via charged-current (CC), 
or neutral-current (NC) interactions, or via both interaction types. The classification 
can be performed by considering the quantum numbers of the final-state fermion pairs. 
In the following, / and F denote different fermions (/ 7^ F) that are neither electrons 
nor electron neutrinos {f,F 7^ e~,t'e), and their weak-isospin partners are denoted by 
/' and F', respectively. We find the following 11 classes of processes (in parenthesis the 
corresponding classification of Ref. [ ^ is given): 

(i) CC reactions: 

(a) e+e- ff'FF', {CCll family), 

(b) e+e- ^ i/ee+//', {CC20 family), 

(c) e+e- ^ ffe-Ue, {CC20 iamily), 
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a) abelian graph 





MPb, CTb) 



fc{Pc, (Tc) 
fdiPd, 

fe{Pe,(Te) 
ffiPf^^f) 



b) non-abelian graph 



fa{Pa, CTa 



fb{Pb, 




ffiPf^^f) 



Figure 1: Generic diagrams for e"'"e —>■ 4/ 



fii) NC reactions: 



(a) e+e- ^ ffFF, 

(b) e+e- ^ ////, 

(c) e+e" e'e^ff, 

(d) e+e" — i> e~e"'"e~e"'", 

(iii) Mixed CC/NC reactions: 

(a) e+e- ffff, 

(b) e+e" Uei^eff, 

(c) e+e" Ue^e^e^e, 

(d) e+e~ — i> i/et'ee~e+, 



(iVC5^ family), 
(iVC^-i 5 family), 
(iVC^5 family), 
(iVC^-.?^ family), 

{mix43 family), 
{NC21 family), 
(iVC^- 5 family), 
{mix56 family). 



The radiation of an additional photon does not change this classification. 

2.3 Generic diagrams and amplitudes for e+e~ — > 4/ 

In order to explain and to illustrate our generic approach we first list the results 
for e+e~ 4/. All these processes can be composed from only two generic diagrams, 
the abelian and non-abelian diagrams shown in Figure All external fermions fa,...j 
are assumed to be incoming, and the momenta and helicities are denoted by Pa,...,f and 
(^a,...ji respectively. The helicity amplitudes of these diagrams are calculated within the 
Weyl-van der Waerden (WvdW) formalism following the conventions of Ref. [ |18[ (see 
also references therein). 

2.3.1 Leptonic and semi-leptonic final states 



We first treat purely leptonic and semi-leptonic final states. In this case, none of the 
gauge bosons in the generic graphs of Figure IT] can be a gluon, and the colour structure 
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trivially leads to a global factor Nj^Nj^, which is equal to 1 or 3, after summing the squared 
amplitude over the colour degrees of freedom. The results for the generic amplitudes are 



MviV2 {Pa:Pb:Pc:Pd:Pe:Pf) 

X ^i(^±p)M?!^±M^;.'...(p„,p,,p„,p,,p.,p,), (2.8) 

J^YWW \Pa. Pb . Pc . Pd -Pe-Pf) 

X Pv{Pa+ Ph)Pw{Pc + Pd)Pw{Pe+ Pf) Al''{.Pa,Ph,Pc,Pd,Pe,Pf), (2-9) 

where the vector-boson propagators are abbreviated by 

Py^P)- 2 \^2 ^ V = A,Z,W,g, Ma = M, = Q. (2.10) 

— My 

(The case of the gluon is included for later convenience.) The auxiliary functions 
and ^43° are expressed in terms of WvdW spinor products, 

At^^{Pa,Pb,Pc,Pd,Pe,Pf) = {PaPc){PbPf)*{{PbPd)*{PbPe) + ipdP f)* ipeP f)) , 

At^'{Pa,Pb,Pc,Pd,Pe,Pf) = A+^^{pa,Pb,Pc,Pd,Pf,Pe), 

At~^{Pa,Pb,Pc,Pd,Pe,Pf) = At^^{pa,Pb,Pd,Pc,Pe,Pf), 

At~~iPa,Pb,Pc,Pd,Pe,Pf) = At^^{pa,Pb,Pd,Pc,Pf,Pe), 

'{Pa,Pb,Pc,Pd,Pe,Pf) = {M~''"~''''{Pa,Pb,Pc,Pd,Pe,Pf)) , (2.11) 



At(Pa,Pb,Pc,Pd,Pe,Pf) = {PbPd)*{PbPf)*{PaPb){PcPe) + {PbPd)* (PdPf)* {PaPe) {PcPd) 

+ {PbPfy{PdPf)*{PaPc){PePf), 
A3(Pa,Pb,Pc,Pd,Pe,Pf) ^ A^(Pb,Pa,Pc,Pd,Pe,Pf)- (2.12) 

The spinor products are defined by 



e-i<^. cos ^ sin ^ - e"''^'' cos ^ sin ^ 
2 2 2 2 



(2.13) 



where pa, Qa are the associated momentum spinors for the momenta 

p'^ — po(l, sin 9p cos (pp, sin 9p sin (f)p, cos 9p) , 

(f — Poi^: sin 6q cos (f)q, sin 9q sin (f)q, cos 6q). (2-14) 

Incoming fermions are turned into outgoing ones by crossing, which is performed by 
inverting the corresponding fermion momenta and helicities. If the generic functions are 
called with negative momenta —p, —q, it is understood that only the complex conjugate 
spinor products get the corresponding sign change. We illustrate this by simple examples: 

^(p, q) = (pq) = ^(p, -q) = ^(-p, q) = ^(-p, -?), 

B{p, q) = {pq)* = -B{p, -q) = -B{-p, q) = B{-p, -q). (2.15) 
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We have checked the resuhs for the generic diagrams against those of Ref. [ 16 1 and found 
agreement. 

Using the results for the generic diagrams of Figure |1], the helicity amplitudes for all 
possible processes involving six external fermions can be built up. It is convenient to 



construct first the amplitudes for the process types CC(a) and NC(a) (see Section |2^ ) 
in terms of the generic functions (|2.8|) and (p^) , because these amplitudes are the basic 
subamphtudes of the other channels. The full amplitude for each process type can be 
built up from those subamphtudes by appropriate substitutions and linear combinations. 
We first list the helicity amplitudes for the CC processes: 



0-+ , 0-_ , Cri ,(T2 , (T3 , 0-4 



'CCa 



Mww ' '''' ''*{p+,P-,-h,-k2,-ks,-k4) 



+ 



,^0"l,^0"2i^o-3,— (T4 



V=-y,Z 



I A cri, — ''"2,o-+,o-_, — (73,— 0-4 
I \A-<^3,-<^4,o-+,0'-,-<Ji,-cr2 

I A ^-'T3,-a-4,— o-l,— cr2,(T+,(7_ 



A ^ 0"+ ,0"_ ,(7l ,(T2 ,0-3 ,0-4 

-'^'CCb 



M 



{p+,p^,ki,k2,k3,k4) 

o-+,cr_ ,0-1,0-2 ,0-3,0-4/ T 

CCa yP+iP--,!^ 



-^CCa""'"""""'"'"*( 



CCc 



{p+,p^,ki,k2,k3,k4 

CCa ''"'''*{p+,p-,ki,k2,k3,k4 

-0-3, 
CCa 



0-|_ ,o-_ ,01 ,0-2 ,03 ,04 

—0-3,0-- ,CT1 ,02 , — CT+ ,04 



The ones for the NC processes are given by 

{p+,P-,h, 



01,— 02,— 03,— 04 



yi,y2=7,^" 



-01,- 



03, — 04,— 01,~02 
— 01, — CT2,CT+,0_ , — 03,-04 



A ^ ^03, — 04,0+,0_ , — 01, — 02 



'1 V2 

-01,— 02, — 03, — 04,0+, O- 



V1V2 



I A/f~''"3!~''"4i~''"li~''"2,0+,0- 



(p+,P_, -fci, -fc2 


7 7 \ 


(-/ei, -k2,p+,p- 


-A;3, -ki) 


( ICa T) , T) 

\ i^3i "'4)f+)f— 




{-ki, -k2, -fcs, - 


h,P+,P~) 


(-^3, -^4, -fcl, - 


h,P+,P-) 


1 1 \ 




, /Cs, /C4), 




/C3, k^) 




,h,-p+,k4). 




(p+,p_, -/ci, -A;2 


-k-i, -ki) 


{p+,P-, -A;3, -A;4 


-ki, -k2) 


{-ki,-k2,p+,p- 


-ks, -ki) 


(-/C3, -k4,p+,p- 


-ki, -^2) 


{-ki, -k2, -ks, - 


h,P+,P-) 


(-^3, -A;4, -^1, - 


k2,p+,p-) 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



0+ ,o_ ,01 ,02 ,03 ,04 



'NCb 



{p+,p-,ki,k2,k3,ki) 

0+ ,o_ ,01 ,02,03,04 , 



-MN"c:-"''^"'^"'^'(p+,P-,fci,fc2,A;3,fc4) 

_^0+,0_,03,02,01,04^^^^^_^^^^^^^^^^^^^ 



(2.20) 
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- Mnc/-'"'^^''"'"'^'(-^i,P-, h, h, h), (2.21) 

= -MN^c'r-'"^'"^''^^''^^(p+,p-, h, h, h, h) 

- -Mnc:-''^^'^^'^^''^^(p+,p-, h, h, h, h) 

- -Mnc/-'~'+''"'"'^'(-^i,P-, ^2, A;3, h) 

+ M^cT-'"'^'''"''"''\-k3,P-, k„ k„ k,). (2.22) 
Finally, the helicity amplitudes for reactions of mixed CC/NC type read 

A^CC/NCa k,, ks, k,) 

= A^n"c':-"^^''^^''^^'''(p+,P-, ku k2, ks, k,) 

- Ml^d:-'''"''"''"'''iP+,P-, ki, k,, ks, k,), (2.23) 



■^CC/N6b "^"''"'^'(P+^P-, ^1, ^2, ^3, ^4) 

= A4^+c:-''^^''^^''^-'^^(p+,p_, k„ k„ ks, k,) 

- MccT''"''''"'-'"'^''''i-k3, -k,, k,, k2), (2.24) 

■MSc/n6?'"''^""(p+,P-, k2, ks, k,) 

= -MNt^':-''^''^'''^^''^'(p+,P-, A:2, ^3, A;4) 

- -Mn"c:-''^^'^^''^^''^^(p+,P-, A;3, A;i, A;4) 

- Mcc'r''-"^^"''"''-"- {-k„ -k„ ks, -P-, k,) 

+ A<cSr''"~''^'''''"~"-(-fcl, -A:4, ^3, -P-, -P+, k2) 
+ -McCa'"'^""'^^''"'^"~'^-(-^3, -k2, k,, -p., -p+, k,) 

- MccT""'^'''"''"~^^i-k3, -k„ k„ -p., k2), (2.25) 



■^St^/N6d ''"'^"''(P+.P-, ^1, k2, ks, k,) 

= MNt^':-''"'^"'^"'^(p+,P-, ^1, ^2, ^3, ^4) 

- MNSr-''^^''^''"'^^''^'(-^3,P-, ^1, A;2, k,) 

- A^ct^'r-''"'^"'^"'^^(P+,P-, ^1, A:4, ^3, ^2) 

+ MSt^r"'""'^-''^"'^'(P+, -^4, k„ A;3, ^2) 

- A4cSr'^"'^""'^-'"'^"''^^(-^3, -^4, A;^). (2.26) 

The relative signs between contributions of the basic subamplitudes A^cca and A^NCa to 
the full matrix elements account for the sign changes resulting from interchanging external 
fermion lines. 
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For the CC reactions, the amphtudes A^cca are the smallest gauge- invariant subset of 
diagrams [ In the case of NC reactions, the amplitudes A^NCa are composed of three 



separately gauge-invariant subamplitudes consisting of the first two lines, the two lines in 
the middle, and the last two lines of ( p.l9| ). 

2.3.2 Hadronic final states 

Next we inspect purely hadronic final states, i.e. the cases where all final-state fermions 
fi are quarks. This concerns only the channels CC(a), NC(a), NC(b), and CC/NC(a) 
given in Section [2.2|. The colour structure of the quarks leads to two kinds of modifications. 



Firstly, the summation of the squared amplitudes over the colour degrees of freedom can 
become non-trivial, and secondly, the possibility of virtual-gluon exchange between the 
quarks has to be taken into account. More precisely, there are diagrams of type (a) in 
Figure |I| in which one of the gauge bosons Vi,2 is a gluon. The other gauge boson of 
Vi^2 can only be a photon or Z boson, since this boson has to couple to the incoming 
e"'"e~ pair. Consequently, there is an impact of gluon-exchange diagrams only for the 
channels NC(a), NC(b), and CC/NC(a), but not for CC(a). This can be easily seen by 
inspecting the generic diagrams in Figure |l|: the presence of a gluon exchange requires 
two quark-ant iquark pairs qq in the final state. 

We first inspect the colour structure of the purely electroweak diagrams. Since the 
colour structure of each diagram contributing to the basic channels CC(a) and NC(a) 
is the same, the corresponding amplitudes factorize into a simple colour part and the 
"colour- singlet amplitudes" A^cca and AlNCa, given in (|2.16| ) and (|2.19| ), respectively. 



The amplitudes for NC(b) and CC/NC(a) are composed from the ones of CC(a) and 
NC(a) in a way that is analogous to the singlet case, but now the colour indices q of the 
quarks fi have to be taken into account. Indicating the electroweak amplitudes for fully 
hadronic final states by "had, ew" , and writing colour indices explicitly, we get 

■^Sca,h;d!;w;?;2*c3,e4(p+.p-' ^2, h, h) 

= -MScr''''''''"''"''(P+,P-, h, k2, h, h) 4,,,5,3,„ (2.27) 

= -MNt^':-''^^''^^"^^''^^(p+,P_, h, h, h) 5eie25c3C4, (2.28) 
■^fKh.had'ew'cuco^.c-i.CAiP+yP-^ ^1' ^2, ^4) 



NCb,had,ew,ci,C2,C3,C4 

^a+,a_,ai,a2,<73,<74j^^_^^^__^ fcl, h, h) Sc.c^S, 



C3C4 



- -MNt:'r-"^"'^"'^"'^^(p+,P-, h, k2, h, h) 5e3e25ciC4, (2-29) 
^^Cc7N6a'Cd™,C2,C3,C4b+' ^1' ^2, ^3, h) 

- Aiccr''^"""'^"'^^(p+,P-, k,, k,, ks, k2) 4iC443C2- (2.30) 



In the calculation of the gluon-exchange diagrams we can also make use of the "colour- 
singlet" result ( p.8|) for the generic diagram (a) of Figure 0, after splitting off the colour 
structure appropriately. Since each of these diagrams involves exactly one internal gluon. 
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exchanged by the two quark hnes, the corresponding matrix elements can be deduced in 
a simple way from the diagrams in which the gluon is replaced by a photon. The gluon- 
exchange contributions to the channels NC(b) and CC/NC(a) can again be composed 
from the ones for NC(a). Making use of the auxiliary function 

_y^a+,a_,a„..,.3,<x4(p^^^_^ A;i, k2, k^) 



g 

-,2 



I A A 

(—(71 

+ ^-^^■-'^^■-'^^'-^^'^+'^-(-^3, -^4, -A;2,p+,p-)l , (2.31) 



+ M7:^^-'''^"^^"--'''^-''\-k^, -k4,p+.p- -ki, -k2) 



where gs = ^/A^Ta^ is the strong gauge coupling, the matrix elements involving gluon 
exchange explicitly read 

•^NCa,hadlgluonX!c2,C3,C4(^'+'-P-' ^1' ^2, /^3, k4) 

= Mg+'"-'"^'"^'"^'"^(p+,P-, h, k2, ks, k,) i K,c,Ksc,, (2.32) 

A/fO-+,'7-.'71><^2,0-3,0-4 / h„ hA 

•'^lNCb,had,gluon,ci,C2,C3,C4VP+)P-) '^l) '^2, «'3, «'4j 

= A^-+--.-i--^.-3.-4(p^^p_^ 1 A»^,^A«3,^ 

- M^+''^-''^^''^^''^^''^^(p+,P-, A;3, k2, k,, k,) \ A»3,,A« (2.33) 

•^CC/NCa'had,gluon,ci,C2,C3,C4(^'+'-P-' ^1' ^2, ^3, ^4) 

= A^g+'"-'"^'"^'"^'"*(p+,P-, A;2, k^, k,) i A^^,^A«3,^. (2.34) 

The colour structure is easily evaluated by making use of the completeness relation 
^ij^ki — ~^^ij^ki + '^^ii^jk for the Gell-Mann matrices A^^-. 

The complete matrix elements for the fully hadronic channels result from the sum of 
the purely electroweak and the gluon-exchange contributions, 

«yO-+,(T_,(Ti,0-2,(73,(T4 <T_|. ,(T_ ,(71 ,(72 ,(73 ,(74 (7+ ,(7_ ,(7l ,(72 ,(73 ,(74 /o q|-\ 

...,had,ci,C2,C3,C4 ...,had,cw,ci,C2,C3,C4 ...,had,gluon,ci,C2,C3,C4" y^.OOj 

The gluon-exchange contributions are separately gauge-invariant. 

For clarity, we explicitly write down the colour-summed squared matrix elements 
for the fully hadronic channels. Abbreviating A4"'+''^-''^°""'^'"''""''^{p+,p-.,ka,kb,kc,kd) by 
Al...(a, b, c, d) we obtain 

E |Mcca,had(l,2,3,4)|2 = 9|Mcca(l,2,3,4)|^ (2.36) 

colour 

E l-MNCa,had(l, 2, 3, 4)|2 = 9|AtNCa(l, 2, 3, 4)1^ + 2|Mg(l, 2, 3, 4)|^ (2.37) 

colour 

E |-MNCb,had(l,2,3,4)|2 



colour 



9|MNCa(l, 2, 3, 4)1^ + 9| A^NCa(3, 2, 1, 4)1^-6 Re {MNCa(l, 2, 3, 4)M;jca(3, 2, 1, 4)} 
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+ 2\M,{1, 2, 3, 4) |2 + 2| A^g(3, 2, 1, 4) p + | Re {M,{1, 2, 3, 4)X*(3, 2, 1, 4)} 
- 8 Re {A<NCa(l, 2, 3, 4)M*(3, 2, 1, 4)} - 8 Re {A<NCa(3, 2, 1, 4)A<*(1, 2, 3, 4)}, 

(2.38) 

E l-^CC/NCa,had(l,2,3,4)|2 
colour 

= 9|A<NCa(l, 2, 3, 4) |2 + 9| A^cca(l, 4, 3, 2)|2 - 6 Re {M^c.{l, 2, 3, i)M*cc.{h 4, 3, 2)} 
+ 2\M^{1, 2, 3, 4) |2 - 8 Re {A^cca(l, 4, 3, 2)A^*(1, 2, 3, 4)}. (2.39) 

Owing to the colour structure of the diagrams, a non-zero interference between purely elec- 
troweak and gluon-exchange contributions is only possible if the four final-state fermions 
can be combined into one single closed fermion line in the squared diagram. This im- 
plies that fermion pairs must couple to different resonances in the electroweak and the 
gluon-exchange diagrams, leading to a global suppression of such interference effects in 
the phase-space integration (see Section ^]5|). 

2.4 Generic functions and amplitudes for e+e" — > 4:f^ 

The generic functions for e"'"e~ —>■ 4/7 can be constructed in a similar way. The idea 
is to combine the contributions of all those graphs to one generic function that reduce 
to the same graph after removing the radiated photon. These combined contributions 
to e^e~ — > 4/7 are classified in the same way as the diagrams for the corresponding 
process e"'"e~ —>■ 4/, i.e. the graphs of Figure |I| also represent the generic functions for 
e"'"e~ —>■ 4/7. Finally, all amplitudes for e"'"e~ —>■ 4/7 can again be constructed from 
only two generic functions. Note that the number of individual Feynman diagrams ranges 
between 14 and 1008 for the various processes. We note that the generic functions can in 
fact be used to construct the amplitudes for all processes involving exactly six external 
fermions and one external photon, such as e~e~ Af'j and e7 — >■ 5/. 

As a virtue of this approach, the so-defined generic functions fulfill the QED Ward 
identity for the external photon, i.e. replacing the photon polarization vector by the 
photon momentum yields zero for each generic function. This is simply a consequence 
of electromagnetic charge conservation. Consequently, in the actual calculation in the 
WvdW formalism the gauge spinor of the photon drops out in each contribution separately. 

Assuming the external fermions as incoming and the photon as outgoing, the generic 
functions read 

j^aa,at,ac,aa,<7e,af,X^Q^^ Qfe, Qc, Qd, Qe, Qf, Pa, Pb, Pc, Pd, Pe,Pf, k) 
-tv c;^„ _^^c;^^ ^Vifafg^Vifgh^Vifcfa^Vifeff 

X A2'''''"'''"'^{Qa,Qb,Qc,Qd,Qe,Qf,Pa,Pb,Pc,Pd,Pe,Pf,k), (2.40) 

Mvww''''''''''"'''^'^{Qa, Qb, Qc, Qd, Qe, Qf, Pa, Pb, Pc, Pd, Pe, Pf, k) 

= -AV2e^6a^-a^Sa„+S^^-6a^,^+6a^- {Qc - Qdjgvwwgvf^f^^dwuf^Swj.ff 

X A^'"^{Qa,Qb,Qc,Qd,Qc,Qf,Pa,Pb,Pc,Pd,Pe,Pf,k), (2.41) 
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with the auxihary functions 

^t^^^{Qa,Qb,Qc,Qd,Qe,Qf,Pa,Pb,Pc,Pd,Pe,Pf,k) = -(PaPc)| 

Pvi (Pc + Pd)Pv2 (Pe + Pf) 

'{PbPfY ( Qc-Qd (PaPc) 



-{{PbPdY{PhPe) + {PdPf)*{PePf)) 



+ 
+ 



(PaA:) \{.Pb+Pe+PfY (Pck) 

+ ^^^-^T^i(Po.Pd)* iPaPe) + iPcPd)* iPcPe))] 

[Pa+Pc+Pdr \Pek) ) 

Qb{{PaPdY(j)aPe) + {P cP d)* {p cP e)) {{PbP f)* {p aPb) " {Pfk)*(j)ak)) 

{Pa+Pc+PdY{Pak){pbk) 
{Qa + Qc- Qd){PbPf)*{PcPd)*{PaPc){{Pbky*{pbPe) " {Pfk)*{PePf)) 



X 



(Pa +Pc+ PdYiPb +Pe+ PfY{Pak) 

Qd-{Qc-Qd)'^{Pd-k)PvAPc+Pd) r> / , /AO r , V \* 

[Pb+Pe+PfY 

{PcPd){{PbPdY{PbPe) + {PdPfY{PePfY + {Pck){{pbkY{PbPe) - {PfkYjPePf)) 

{Pck){pdk) 

, Qf - {Qe - Qf)2iPf ■ k)Pv,{pe + Pf) „ . ^ . ^ ,x 

+ Pv^iPc+P^)Pv.iPe+Pf - k) 

^ {{PbPfYjPePf) + {PbkY{Pek))i{PaPdY{PaPe) + {PcPdY jPcPe)) 

(pek) (pfk) 

Pa,Pb,Pc,Pd,Pe,Pf,k) 
= Qft, Qc, Qd, -Qf: -Qe:Pa:Pb,Pc,Pd:Pf:Pe: k), 

^^i^Qa-iQb-iQc-iQd-iQe-iQf-i Pa,Pb,Pc,Pd,Pe,Pf,k) 

— ^2'^'^'^ {Qai Qbi ~Qdi ~Qci Qei Q f 1 Pa,Pb,Pd,Pc,Pe,Pf,k), 
■^2 '^i,Qa-i Qbi Qci Qdi Qei Qfi PaiPbiPciPdiPeiPfik) 

= >lj+++((5a, Qb, -Qdi -Qci -Qfi -QeiPaiPbiPdiPciPfiPei k), 
■^2 {Qai Qbi Qci Qdi Qei Qfi 

PaiPbiPciPdiPeiPfik) 
= At^'^'^iQbi Qai -Qei -Qfi -Qci -QdiPbiPaiPeiPfiPciPdi k), 
■^2 {Qai Qbi Qci Qdi Qei Qfi 

PaiPbiPciPdiPeiPfik) 

— ^2'^'^'^ {Qbi Qai Qfi Qei Qci Qdi PbiPaiPfiPeiPciPdik), 
■^2 {Qai Qbi Qci Qdi Qei Qfi 

PaiPbiPciPdiPeiPfik) 
= Qa, -Qei -Qfi Qdi QciPbiPaiPeiPfiPdiPci k), 

■^2 {Qai Qbi Qci Qdi Qei Qfi PaiPbiPciPdiPeiPfik) 

= ^t~^~^~^{Qbi Qai Qfi Qei Qdi Qci PbiPaiPfiPeiPdiPcik), 

^^-<^-^<^'-(g„, Qf,, g„ Qa, Q„ Qf,PaiPbiPciPdiPeiPfi k) (2.42) 

= [A^"" {QaiQbiQciQdiQeiQfiPaiPbiPciPdiPeiPfik)) 
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and 



^t^iQa, Qb, Qc, Qd, Qe,Qf, Pa, Pb, Pc, Pd, Pe,Pf, k) 

= PviPa + Pb)PwiPc + Pd)PwiPe + ^/ ) 

X ( {PbPdY{PbPfY{PaPb) {PcPe) + {PbPd)*{PdPfY{PaPe) {PcPd) 
+ {PbPfY{PdPfY{PaPc){PePf)) 

+ Pv{,Pa +Pb- k)Pw{Pc + Pd)PwiPe + Pf) (^p k)\pbk) 

X{ {PdPfYi {PaPe){PcPd){{PbPdY{PaPb) - {PdkYiPak)) 
+ {PaPc){PePf){{PbPfY{PaPb) " {PfkYiPakY] 
+ {PcPe){{PbPdY{PaPb) - {PdkY{Pak)){{pbPfY{PaPb) " {pfkY{PakY} 

,p/ , NO / , i\p f , .Qd-{Qc-Qd)2{pd-k)Pw{Pc + Pd) 

+ PvKPa + Pb)Pw[Pc +Pd- k)Pw[Pe + Pf) (p k) (p^k) 

X{ {PbPfYl {PaPc){PePf){{PdPfY{PcPd) - {PfkYiPck)) 
+ {PcPe){PaPb){{PbPdY{PcPd) + {PbkYiPck))] 
+ {PaPe){{PdPfY{PcPd) - {PfkY{Pck)){{pbPdY{PcPd) + {PbkYiPck))} 

X{ {PbPdYl {PcPe){PaPb){{PbPfY{PePf) + {PbkYiPek)) 
+ {PaPe){PcPd){{PdPfY{PePf) + {PdkY {Pek))] 
+ {PaPc){{PbPfY{PePf) + {PbkY{Pek)){{pdPfY{PePf) + {PdkY {Pek))} , 
^S^iQa, Qb, Qc, Qd, Qe, Qf,Pa, Pb, Pc, Pd, Pe,Pf, k) 

= A^+{-Qb, -Qa, Qc, Qd, Qe, Qf,Pb, Pa, Pc Pd, Pe, Pf, k), 
^S"'" (Qa, Qb, Qc, Qd, Qc, Qf, Pa, Pb, Pc, Pd, PcPf, k) (2.43) 
= (A^'''''~^{Qa,Qb,-Qd,-Qc,-Qf,-Qe,Pa,Pb,Pd,Pc,Pf,Pc,k)) r , a* r ■^■ 

The replacements Py — > Py after the complex conjugation in the last lines of (|2.43| ) and 
( p.44| ) ensure that the vector-boson propagators remain unaffected. Note that the vector- 
boson masses do not enter explicitly in the above results, but only via Py. In gauges 
such as the 't Hooft-Feynman or the unitary gauge this feature is obtained only after 
combining different Feynman graphs for e"'"e~ 4:f'y', in the non-linear gauge (p.3|) this 
is the case diagram by diagram. 

The helicity amplitudes for e"'"e~ — > 4/7 follow from the generic functions A^viVi and 
M-vww of ( |2.41| ) in exactly the same way as described in Section ^]B| for e"'"e~ — > 4/. This 
holds also for the gluon-exchange matrix elements and for the colour factors. Moreover, 
the classification of gauge-invariant sets of diagrams for e"^e~ —>■ 4/ immediately yields 
such sets for e"'"e~ — >■ 4/7, if the additional photon is attached to all graphs of a set in all 
possible ways. 

We have checked analytically that the electromagnetic Ward identity for the external 
photon is fulfilled for each generic contribution separately. In addition, we have numeri- 
cally compared the amplitudes for all processes with amplitudes generated by Madgraph 
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[ 25 1 for zero width of the vector bosons and found complete agreement. We could not 
compare our results with Madgraph for finite width, because Madgraph uses the unitary 
gauge for massive vector-boson propagators and the 't Hooft-Feynman gauge for the pho- 
ton propagators, while we are using the non-linear gauge (|2.3| ). Therefore, the matrix 
elements differ after introduction of finite vector-boson widths. While the calculation 
with Madgraph is fully automized, in our calculation we have full control over the matrix 
element and can, in particular, investigate various implementations of the finite width. 

A comparison of our results with those of Refs. [ |10|, |ll| , which include only the matrix 
elements that involve two resonant W bosons, immediately reveals the virtues of our 
generic approach. 

2.5 Implementation of finite gauge-boson widtiis 

We have implemented the finite widths of the W and Z bosons in different ways: 

• fixed width in all propagators: Pvijp) = [p"^ ~ + i^yry]"^, 

• running width in time-like propagators: Pv{p) = [p'^ ~ My + ip^(rv/Mv)^^(j5^)]~"'^, 



• complex-mass scheme: complex gauge-boson masses everywhere, i.e. y My — iMyFy 
instead of My in the propagators and in the couplings. This results, in particular, 
in a constant width in all propagators, 

Pv{p) = [p' -M^ + iMyTv]-\ (2.44) 
and in a complex weak mixing angle: 

" " " M| - iMzFz ■ ^ ' 

The virtues and drawbacks of the first two schemes have been discussed in Ref . [ ^ . 
Both violate SU(2) gauge invariance, the running width also U(l) gauge invariance. The 
complex-mass scheme obeys all Ward identities^ and thus gives a consistent description 
of the finite-width effects in any tree-level calculation. While the complex-mass scheme 
works in general, it is particularly simple for e"'"e~ 4/7 in the non-linear gauge ( p.3|) . 
In this case, no couplings involving explicit gauge-boson masses appear, and it is sufficient 
to introduce the finite gauge-boson widths in the propagators [cf. (|2.44| )] and to introduce 



the complex weak mixing angle (|2.45| ) in the couplings. We note that a generalization of 



this scheme to higher orders requires to introduce complex mass counterterms in order 
to compensate for the complex masses in the propagators [|2^. We did not consider the 
fermion-loop scheme [ which is also fully consistent for lowest-order predictions, 

since it requires the calculation of fermionic one-loop corrections to e"'"e~ — > 4/7 which is 
beyond the scope of this work. 



^In the context of electromagnetic gauge invariance, the introduction of a complex mass has been 
proposed in Ref. [ M 
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3 The Monte Carlo programs 



The cross section for e^e 
(27r)4-3" 



da 



2s 



— ^ 4/(7) is given by 

n 



P++P- 



i=l 



X \M{p+,p-, ki,...,k„ 



(3.1) 



where n = 4, 5 is the number of outgoing particles. The helicity amplitudes Ai for 
e"'"e~ — > 4/(7) have been calculated in Sections and |2]J, respectively. The phase- 
space integration is performed with the help of a Monte Carlo technique, since the Monte 
Carlo method allows us to calculate a variety of observables simultaneously and to easily 
implement cuts in order to account for the experimental situation. 

The helicity amplitudes in (|3.1| ) exhibit a complicated peaking behaviour in different 
regions of the integration domain. In order to obtain a numerically stable result and to 
reduce the Monte Carlo integration error we use a multi-channel Monte Carlo method [ 
TBI, Ol, which is briefly outlined in the following. 

Before turning to the multi-channel method, we consider the treatment of a single 
channel. We choose a suitable set $ of 3n — 4 phase-space variables to describe a point in 
phase space, and determine the corresponding physical region V and the relation fci($) 
between the phase-space variables $ and the momenta fci, . . . The phase-space inte- 
gration of (|3.1|) reads 



da = J^d$p(k,{$))f[ki{$ 
(27r)4-3n 



(3.2) 



2s 



M {p+,p.,ki{$),...,kn{$) 



where p is the phase-space density. For the random generation of the events, we further 
transform the integration variables $ to 3n — 4 new variables r = (fj) with a hypercube 
as integration domain: $ = h{f) with < < 1. We obtain 



V 



d$p(fc,($))/(fc,($) 



1 /(fc,(Mf)) 

Civ — ^ 

g{ki{h{f))) 



(3.3) 



where g is the probability density of events generated in phase space, defined by 



1 



dh{f 



df 



(3.4) 



r=/i-i($) 



If / varies strongly, the efficiency of the Monte Carlo method can be considerably enhanced 
by choosing the mapping of random numbers r into $ in such a way that the resulting 
density g mimics the behaviour of |/|. For this importance sampling, the choice of $ 
is guided by the peaking structure of /, which is determined by the propagators in a 
characteristic Feynman diagram. 

We choose the variables $ in such a way that the invariants corresponding to the prop- 
agators are included. Accordingly, we decompose the n-particle final state into 2^2 
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scattering processes with subsequent 1^2 decays. The variables $ consist of Lorentz 
invariants Si,ti, defined as the squares of time- and space-hke momenta, respectively, and 
of polar and azimuthal angles 6i,4>i, defined in appropriate frames. A detailed description 
of the parameterization of an n-particle phase space in terms of invariants and angles 
can be found in Ref. [ 0. The parameterization of the invariants Si,ti in $ = h{r) is 
chosen in such a way that the propagator structure of the function / is compensated by 
a similar behaviour in the density g. More precisely, if / contains Breit-Wigner reso- 
nances or distributions like s^'^, which are relevant for massless propagators, appropriate 
parameterizations of Sj are given by: 

• Breit-Wigner resonances: 

Si = + MyTy tan[|/i + (t/a - yl)r^] (3.5) 
With ^1,2 = arctan I ^^^^ 1 ; 



1 1 i/Ci-i') 



• propagators of massless particles: 

T^T^l- Si= [s'^^^r, + s'^i^il - ri) 
u = 1 : Si = exp [ln(smax)?^i + ln(smin)(l - ri)] . (3.6) 

The parameter u can be tuned to optimize the Monte Carlo integration and should be 
chosen > 1. The naive expectation u = 2 is not necessarily the best choice, because the 
propagator poles of the differential cross section are partly cancelled in the collinear limit. 
The remaining variables in $ = h{f), i.e. those for which / is expected not to exhibit a 
peaking behaviour, are generated as follows: 

Si = Smax^^i + Smm(l - ^^i), 0i = 27rri, COS 6'i = 2ri - 1. (3.7) 

The absolute values of the invariants ti are generated in the same way as Si. The result- 
ing density g of events in phase space is obtained as the product of the corresponding 
Jacobians, as given in (|3.4| ). In the appendix, we provide an explicit example for an event 
generation with a specific choice of mappings /ci($) and h{f), and for the calculation of 
the corresponding density g. 

The differential cross sections of the processes e+e" 4/ and especially e"'"e~ — > 
4/7 possess very complex peaking structures so that the peaks in the integrand /($) in 
(p.3|) cannot be described properly by only one single density (?($). The multi-channel 



approach [ |T6|, ^ suggests a solution to this problem. For each peaking structure we 
choose a suitable set and accordingly a mapping of random numbers rj into ^k'- 
= hkif) with < Tj < 1, so that the resulting density gu describes this particular 
peaking behaviour of /. All densities gt are combined into one density (?tot that is expected 
to smoothen the integrand over the whole phase-space integration region. The phase-space 
integral of (|3.3|) reads 



'--•^^ ^ ^ ^ ^ gtot[ki{^k)) fc=i-^o gtot[h{hk{r))) 



k=l ' 
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with 



M 



1=1 



{hi^k)) 



9i 



dhiif 



df 



. (3.9) 



r=h-\^k) 



The different mappings hk{f) are called channels, and M is the number of all channels. 

In order to reduce the Monte Carlo error further, we adopt the method of weight 
optimization of Ref. [ ^ and introduce a-priori weights ak,k = 1, . . . , M (a^ > and 
Sfcli Oik = !)• The channel k that is used to generate the event is picked randomly with 
probability ak, i.e. 



In = Y.(^k d$k Pk{ki{^k))gk{ki{^k] 



k=l 
"1 



fi'tot 



M 



dro J2^i^o-/3k-^i)0i/3k 



k=l 
M 



V 



d<^kpk{ki{<^k))gk{ki{<^ 



rl ™ rl 

/ dro J2 0(^0- /3k-i)0i/3k-ro) df 
Jo Jo 



f{ki{hk{r))) 



fjhi^k)) 

gtot{h{^k)) 



gtot{ki(hk{f))y 

where /5o = 0, 13 j = Y.k=i cxk, j = I, ■ ■ ■ , M - 1, Pm = Efcli «fc = 1, and 



(3.10) 



M 



gtotihi^^k)) =J2^i9iikii^ 



(3.11) 



1=1 



is the total density of the event. 

For the processes e"'"e~ —>■ 4/ we have between 6 and 128 different channels, for e^e~ — > 
4/7 between 14 and 928 channels. Each channel smoothens a particular combination of 
propagators that results from a characteristic Feynman diagram. We have written phase- 
space generators in a generic way for several classes of channels determined by the chosen 
set of invariants The channels within one class differ in the choice of the mappings 

(^.5[), (|3.6D, and ( |3.7D and the order of the external particles. We did not include special 
channels for interference contributions. 

The afc- dependence of the quantity 



W{a) 



1 ^ 

— ^[w;(r^,f^)]^ 



(3.12) 



where w = f / gtot is the weight assigned to the Monte Carlo point (rg, f^) of the jth event, 
can be exploited to minimize the expected Monte Carlo error 



5L 



with the Monte Carlo estimate of 



N 



1 ^ 

— ^w(r^,f^) 



(3.13) 



(3.14) 
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by trying to choose an optimal set of a-priori weights. We perform the search for an 
optimal set of by using an adaptive optimization method, as described in Ref. [ . 



After a certain number of generated events a new set of a-priori weights af" is calculated 
according to 



1 A gk{Hhk{r'))) [wirir:>W 



, V ^ ; J : Tar = l. (3.15) 

Based on the above approach, we have written two independent Monte Carlo programs. 
While the general strategy is similar, the programs differ in the explicit phase-space 
generation. 

4 Numerical results 

If not stated otherwise we use the complex-mass scheme and the following parameters: 

a = 1/128.89, as = 0.12, 

Mw = 80.26 GeV, Tw = 2.05 GeV, 

Mz = 91.1884 GeV, rz = 2.46GeV. (4.1) 



In the complex-mass scheme, the weak mixing angle is defined in ( |2.45| ), in all other 
schemes it is fixed by Cw = Mw/Mz, = 1 — c^. 

The energy in the centre-of-mass (CM) system of the incoming electron and positron 
is denoted by ^/s. Concerning the phase-space integration, we apply the canonical cuts 
of the ADLO/TH detector, 

^(Z, beam) > 10°, 9{IJ') > 5°, 9{l,q)>5°, 

^(7,beam) > 1°, ^(7,/)>5°, ^(7,g)>5°, 

E^>0.1 GeV, El > 1 GeV, Eg > 3 GeV, 

m(g, g') > 5 GeV, (4.2) 

where 9{i,j) specifies the angle between the particles i and j in the CM system, and /, 
q, 7, and "beam" denote charged leptons, quarks, photons, and the beam electrons or 
positrons, respectively. The invariant mass of a quark pair qq' is denoted by m(g, q'). The 
cuts coincide with those defined in Ref. [ ^ , except for the additional angular cut between 



charged leptons. The canonical cuts exclude all coUinear and infrared singularities from 
phase space for all processes. 

Although our helicity amplitudes and Monte Carlo programs allow for a treatment 
of arbitrary polarization configurations, we consider only unpolarized quantities in this 
paper. 

All results are produced with 10^ events. The calculation of the cross section for 
e"'"e~ — > e^e~/i"'"/i~ requires about 50 minutes on a DEC ALPHA workstation with 500 
MHz, the calculation of the cross section for e+e~ e"'"e~/i+/i~7 takes about 5 hours. 
The results of our two Monte Carlo programs agree very well. The numbers in parentheses 
in the following tables correspond to the statistical errors of the results of the Monte Carlo 
integrations. 
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4.1 Comparison with existing results 



In order to compare our results for e"'"e~ — ^ 4/ with Tables 6-8 of Ref. [ ^| 
the corresponding set of phase-space cuts and input parameters, i.e. the canonical cuts 
defined in ( [4.2| ), a CM energy of ^/s = 190 GeV, and the parameters a = a(2Mw) = 
1/128.07, as = 0.12, Mw = 80.23 GeV, Tw = 2.0337 GeV, Mz = 91.1888 GeV, and 
Tz = 2.4974 GeV. The value of s^, which enters the couplings, is calculated from 
a(2Mw)/(24) = G^M^/{n^/2) with = 1.16639 x lO-^GeV'^. 

In Table |1], we list the integrated cross sections for various processes e"'"e~ — 4/ with 
running widths and constant widths, and for the corresponding processes e"'"e~ 4/7 
with constant widths. For processes involving gluon-exchange diagrams we give the cross 
sections resulting from the purely electroweak diagrams and those including the gluon- 
exchange contributions. The latter results include also the interference terms between 
purely electroweak and gluon-exchange diagrams. In Table |l| we provide a complete list 
of processes for vanishing fermion masses. All processes e"^e~ — > 4/(7) not explicitly 
listed are equivalent to one of the given processes. 

For NC processes e+e~ — > 4/ with four neutrinos or four quarks in the final state 
we find small deviations of roughly 0.2% between the results with constant and running 



widths. Assuming that a running width has been used in Ref. [|3T|, we find very good 
agreement. 



Unfortunately we cannot compare with most of the publications [ |T^, |T^, [Tj, |T5 
for the bremsstrahlung processes e"'"e~ — >■ 4/7. In those papers, either the cuts are not 
(completely) specified, or collinear photon emission is not excluded, and the corresponding 
fermion-mass effects are taken into account. Note that the contributions of collinear 
photons dominate the results given there. 



We have compared our results with the ones given in Refs. [|T^, |TI[, where the total 
cross sections for e"''e~ Af'j have been calculated for the purely leptonic and the semi- 
leptonic final states. As done in Refs. [ |l^, [ill] only diagrams involving two resonant 



W bosons have been taken into account for this comparison. Table ^ contains our results 
corresponding to Table || of Ref. [ |10|- Based on Refs. [ 0, |Tl[], we have chosen a/s = 



200 GeV and the input parameters a = 1/137.03599, Mw = 80.9 GeV, Fw = 2.14 GeV, 
Mz = 91.16 GeV, Fz = 2.46 GeV, obtained from a/{2sl) = G^,M^/{Tiy/2) with 
G^ = 1.16637 X 10^^ GeV^^, and constant gauge-boson widths. The energy of the photon 
is required to be larger than -E^,min, and the angle between the photon and any charged 
fermion must be larger than 6'^,niin- A maximal photon energy is required, < 60 GeV, 
in order to exclude contributions from the Z resonance. Our results are consistent with 



those of Refs. [ |T0], ^ within the statistical error of 1% given there. In some cases we 



find deviations of 2%.0 

4.2 Comparison of finite-width schemes 

As discussed in Refs. [ particular care has to be taken when implementing the 



finite gauge-boson widths. Differences between results obtained with running or constant 
widths can already be seen in Table |l], where a typical LEP2 energy is considered. In 

^Note that the input specified in Refs. [ |l^, |l^ is not completely clear even if the information of both 
publications is combined. 
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rr / fh 
u 1 LD 


e+e- ^ 4/ 
running width 


e+e- 4/ 
constant width 


e+e — s> 4/7 
constant width 




256.7(3) 


257.1(7) 


89.4(2) 




227.4(1) 


227.5(1) 


79.1(1) 




228.7(1) 


228.8(1) 


81.0(2) 




218.55(9) 


218.57(9) 


76.7(1) 


e~e+e~e+ 


109.1(3) 


109.4(3) 


38.8(4) 


e-e+//-//- 


116.6(3) 


116.1(3) 


13.1(1) 




5.478(5) 


5.478(5) 


3.37(1) 




11.02(1) 


11.02(1) 


6.78(3) 




14.174(9) 


14.150(9) 


5.36(1) 




17.78(6) 


17.73(6) 


6.63(2) 




10.108(8) 


10.103(8) 


4.259(9) 




4.089(1) 


4.082(1) 


0.7278(7) 




8.354(2) 


8.337(2) 


1.512(1) 




4.069(1) 


4.057(1) 


0.7434(7) 




8.241(2) 


8.218(2) 


1.511(1) 


ude~Pe 


693.5(3) 


693.6(3) 


220.8(4) 




666.7(3) 


666.7(3) 


214.5(4) 


e~e+uu 


86.87(9) 


86.82(9) 


32.3(2) 


e~e"'"dd 


43.02(4) 


42.95(4) 


16.17(8) 




24.69(2) 


24.69(2) 


12.70(4) 


d d 


23.73(1) 


23.73(1) 


10.43(2) 


Z/e^'cU u 


24.00(2) 


23.95(2) 


6.84(1) 


Z/e^'ed d 


20.657(8) 


20.62(1) 


4.319(6) 




21.080(5) 


21.050(5) 


6.018(9) 


d d u^u^ 


19.863(5) 


19.817(5) 


4.156(5) 


uudd 


2064.1(9), 2140.8(9) 


2064.3(9), 2141(1) 


615(1), 672(1) 


u d s c 


2015.2(8) 


2015.3(8) 


598(1) 


uuuu 


25.738(7), 71.28(4) 


25.721(7), 71.30(4) 


9.78(2), 42.1(1) 


dddd 


23.494(6), 51.35(3) 


23.448(6), 51.32(3) 


5.527(7), 28.68(4) 


uucc 


51.61(1), 144.72(9) 


51.57(1), 144.75(9) 


19.61(4), 86.1(2) 


uuss 


49.68(1), 126.52(8) 


49.62(1), 126.52(8) 


15.17(2), 75.1(2) 


ddss 


47.13(1), 104.79(6) 


47.02(1), 104.74(6) 


11.10(2), 59.2(1) 



Table 1: Integrated cross sections for all representative processes e+e^ 4/ with running 
widths and constant widths and for the corresponding processes e+e^ — > 4/7 with con- 
stant widths. If two numbers are given, the first results from pure electroweak diagrams 
and the second involves in addition gluon-exchange contributions. 
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IGeV 


5GeV 


lOGeV 


15GeV 




^7, mill 


a/fb 


leptonic 
process 


1° 


53.54(8) 


27.57(3) 


16.96(2) 


11.22(2) 


5° 


32.65(4) 


16.98(3) 


10.48(2) 


6.94(1) 


10° 


23.48(3) 


12.30(2) 


7.61(2) 


5.04(1) 


15° 


18.03(2) 


9.51(2) 


5.90(1) 


3.90(1) 


semi-leptonic 

process 


r 


141.9(2) 


71.90(8) 


43.56(5) 


28.26(4) 


5° 


86.8(1) 


44.25(6) 


26.78(4) 


17.40(3) 


10° 


62.29(7) 


31.92(5) 


19.40(4) 


12.61(3) 


15° 


47.42(5) 


24.50(4) 


14.97(3) 


9.77(2) 



Table 2: Comparison with Table 2 of Ref. [ T^: Cross sections resulting from diagrams 
involving two resonant W bosons for purely leptonic and semi-leptonic final states and 
several photon separation cuts 



(j/fb 


v^ = 


189 GeV 


500 GeV 


2TeV 


lOTeV 


e"'"e~ — > u d /i~z/|j 


constant width 


703.5(3) 


237.4(1) 


13.99(2) 


0.624(3) 


running width 


703.4(3) 


238.9(1) 


34.39(3) 


498.8(1) 


complex-mass scheme 


703.1(3) 


237.3(1) 


13.98(2) 


0.624(3) 


e"'"e~ — > ud/i~i>^7 


constant width 


224.0(4) 


83.4(3) 


6.98(5) 


0.457(6) 


running width 


224.6(4) 


84.2(3) 


19.2(1) 


368(6) 


complex-mass scheme 


223.9(4) 


83.3(3) 


6.98(5) 


0.460(6) 


e+e" — *• ude~z/e 


constant width 


730.2(3) 


395.3(2) 


211.0(2) 


32.38(6) 


running width 


729.8(3) 


396.9(2) 


231.5(2) 


530.2(6) 


complex-mass scheme 


729.8(3) 


395.1(2) 


210.9(2) 


32.37(6) 


e+e" — > ude~z/e7 


constant width 


230.0(4) 


136.5(5) 


84.0(7) 


16.8(5) 


running width 


230.6(4) 


137.3(5) 


95.7(7) 


379(6) 


complex-mass scheme 


229.9(4) 


136.4(5) 


84.1(6) 


16.8(5) 



Table 3: Comparison of different width schemes for several processes and energies 
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y/^l GeV 


189 


500 


V1V2 


WW 


ZZ 


7Z 


77 


WW 


ZZ 


7Z 


77 




26.3 


6.5 


72.5 


94.5 


224 


217 


242 


250 



Table 4: Photon energies ii^y^^^ corresponding to thresholds 

Table |^ we compare predictions for integrated cross sections obtained by using a constant 
width, a running width, or the complex-mass scheme for several energies. We consider 
two semi-leptonic final states for e"''e~ — > 4/(7). The numbers show that the constant 
width and the complex-mass scheme yield the same results within the statistical accuracy 
for e"'"e~ — > 4/ and e"'"e~ — *■ 4/7. In contrast, the results with the running width produce 
totally wrong results for high energies. The difference of the running width with respect 
to the other implementations of the finite width is up to 1% already for 500 GeV. Thus, 
the running width should not be used for linear-collider energies. As already stated above, 
our default treatment of the finite width is the complex-mass scheme in the following. 

4.3 Survey of photon-energy spectra 

In Figure |^ we show the photon-energy spectra of several processes for the typical 
LEP2 energy of 189 GeV and a possible linear-collider energy of 500 GeV. The upper plots 
contain CC and CC/NC processes, the plots in the middle and the lower plots contain 
NC processes. Several spectra show threshold or peaking structures. These structures 
are caused by diagrams in which the photon is emitted from the initial state. The two 
important classes of diagrams are shown in Figure 

The first class, shown in Figure ^a, corresponds to triple-gauge-boson-production sub- 
processes which yield dominant contributions as long as the two virtual gauge bosons 
Vi and V2 can become simultaneously resonant. If the real photon takes the energy E^, 
defined in the CM system, only the energy v^, with 



s' = s-2^E^, (4.3) 

is available for the production of the gauge-boson pair V1V2. If at least one of the gauge 
bosons is massive, and if the photon becomes too hard, the two gauge bosons cannot be 
produced on shell anymore, so that the spectrum falls off for above the corresponding 
threshold E}^^^^. Using the threshold condition for the on-shell production of the V1V2 
pair, 

V^>My, +My„ (4.4) 

the value of E'^^^'^ is determined by 



s - (My, + Mv^ 



\2 



7 



2v^ 



The values of the photon energies that cause such thresholds can be found in Table |^. 
The value i?^'*' corresponds to the upper endpoint of the photon-energy spectrum, which 
is given by the beam energy ^fsjl. Since is fully determined by s and E^^ the 
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Figure 2: Photon-energy spectra for several processes and for = 189 GeV and 500 GeV 
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Figure 3: Diagrams for important subprocesses, where Vi, V2 = W, Z,7, and V3 = 7,g 



contribution of the ViV2-production subprocess to the spectrum quahtatively follows 
the energy dependence of the total cross section for V1V2 production (cf. Ref. [ Fig. 1) 
above the corresponding thresholds. The cross sections for 77 and 7Z production strongly 
increase with decreasing energy, while the ones for ZZ and WW production are comparably 
flat. Thus, the 77 and 7Z-production subprocesses introduce contributions in the photon- 
energy spectra with resonance-like structures, whereas the ones with ZZ or WW pairs yield 
edges. 

The second class of important diagrams, shown in Figure corresponds to the pro- 
duction of a photon and a resonant Z boson that decays into four fermions. These diagrams 
are important if the gauge boson V3 is also resonant, i.e. a photon or a gluon with small 
invariant mass. In this case, the kinematics fixes the energy of the real photon to 



= Ef = (4.6) 



which corresponds to the 7Z threshold in Table ^. This subprocess gives rise to resonance 
structures at -E^^, which are even enhanced by ttg/a in the presence of gluon exchange. 

In the photon-energy spectra of Figure ^ all these threshold and resonance effects are 
visible. The effect of the 7Z peak can be nicely seen in different photon-energy spectra, 
in particular in those where gluon-exchange diagrams contribute (cf. also Figure]^). The 
effect of the WW threshold is present in the upper two plots of Figure ^. In the plot 
for y/s = 189 GeV the threshold for single W production causes the steep drop of the 
spectrum for the pure CC processes above 70 GeV. Note that the CC cross sections are 
an order of magnitude larger than the NC cross sections if the WW channel is open. The 
ZZ threshold is visible in the middle and lower plots for ^/s = 500 GeV. The 7Z threshold 
(resulting from the graphs of Figure ^) is superimposed on the 7Z peak (resulting from 
the graphs of Figure ^b) and therefore best recognizable in those channels where the 7Z 
peak is absent or suppressed, i.e. where a neutrino pair is present in the final state or 
where at least no gluon-exchange diagrams contribute. Processes with four neutrinos in 
the final state do not involve photonic diagrams and are therefore small above the ZZ 
threshold. The effects of the triple-photon-production subprocess appear as a tendency 
of some photon-energy spectra to increase near the maximal value of Ej for two charged 
fermion-antifermion pairs in the final state. 
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189 GeV 500 GeV 2TeV 
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E^l GeV 

Figure 4: Photon-energy spectra resulting from the triple-gauge-boson-production sub- 
processes compared to those resulting from all diagrams (V1V2 includes ZZ, 7Z, and 77) 



4.4 Triple-gauge-boson-production subprocesses 

In Figure ^ we compare predictions that are based on the full set of diagrams with 
those that include only the graphs associated with the triple-gauge-boson-production 
subprocesses, i.e. the graphs in Figure In addition we consider the contributions 
of the ZZ7-production subprocess alone. For CC processes, the photon-energy spectra 
resulting from the W"'"W~7-production subprocess are close to those resulting from all 
diagrams at LEP2 energies, but large differences are found for higher energies and e^ 
in the final state. Note that the spectra are shown on a logarithmic scale. Even at 
LEP2 energies the differences between the predictions for different final states may be 
important, as can be seen, for instance, in Table |l| by comparing the cross sections of 
e"'"e~ — > udyU~i/^7 and e"'"e~ — >• ude'i^eT- In the case of NC processes, already for 
189 GeV the contributions from ZZ7, Z77, and 777 production are not sufficient: in the 
vicinity of the 7Z peak sizeable contributions result from the 7Z-production subprocess 
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a/fb 


ew and gluon 


purely ew 


gluon 


interference 


e"'"e~ — > uucc 


52.98(4) 


21.560(6) 


31.38(3) 


0.04(5) 


e"'"e~ — > UUCC7 


29.8(1) 


10.38(4) 


19.6(1) 


-0.1(1) 


e"'"e~ — i> uuuu 


26.25(2) 


10.765(3) 


15.34(1) 


0.14(2) 


e"'"e~ — > UUUU7 


14.83(7) 


5.16(2) 


9.52(5) 


0.15(9) 


e"'"e~ — i> dduu 


901.2(6) 


876.4(5) 


24.24(2) 


0.6(8) 


e"'"e~ — *■ dduu7 


290(1) 


275(1) 


14.82(8) 


0(1) 



Table 5: Full lowest order cross section (ew and gluon) and contributions of purely elec- 
troweak diagrams (ew), of gluon-exchange diagrams (gluon), and their interference for 
500 GeV 

(Figure Pd) even for the /i"''/i~uu7 final state. For e+e~ e^e+uii7 other diagrams 
become dominating everywhere. The contribution of ZZ7 production is always small 
and could only be enhanced by invariant-mass cuts. Note that the triple-gauge-boson- 
production diagrams form a gauge-invariant subset for NC processes, while this is not the 
case for CC processes. 

4.5 Relevance of gluon-exchange contributions 

In the analytical calculation of the matrix elements for e^e^ —>■ 4/(7) in Section |]we 
have seen that NC processes with four quarks in the final state involve, besides purely 
electroweak, also gluon-exchange diagrams. Table |^ illustrates the impact of these dia- 
grams on the integrated cross sections for a CM energy of 500 GeV. The results for the 
interference are obtained by subtracting the purely electroweak and the gluon contribution 
from the total cross section. For pure NC processes the contributions of gluon-exchange 
diagrams dominate over the purely electroweak graphs. This can be understood from 
the fact that the gluon-exchange diagrams are enhanced by the strong coupling constant. 



and, as discussed in Section ^73| , that the diagrams with gluons replaced by photons 
yield a sizeable contribution to the cross section. For the mixed CC/NC processes the 
purely electroweak diagrams dominate the cross section. Here, the contributions from 
the W"'"W~7-production subprocess are large compared to all other diagrams, even if the 
latter are enhanced by the strong coupling. At 500 GeV the gluon-exchange diagrams 
contribute to the cross section at the level of several per cent. The interference contribu- 
tions are relatively small. As discussed at the end of Section ^.3.2|, this is due to the fact 



that interfering electroweak and gluon-exchange diagrams involve different resonances. 
Note that the interference vanishes for e"'"e~ — > uucc7, and the corresponding numbers 
in Table ^ are only due to the Monte Carlo integration error. 

In Figure |^ we show the photon-energy spectra for the processes e"'"e~ — > u u d d 7 
and e+e" uuuu7 together with the separate contributions from purely electroweak 
and gluon-exchange diagrams. The pure electroweak contributions are similar to the ones 
for e"'"e~ — > ud/i~i/^7 and e"'"e~ — *• e~e"'"uii7 in Figure]^. For the NC process e"'"e~ —>■ 
UUUU7, the photon-energy spectrum is dominated by the gluon-exchange contribution, 
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Figure 5: Electroweak and gluon-exchange contributions to the photon-energy spectra for 
e"'"e~ — > u u d d 7 and e"'"e~ — > u li u u 7 at i/i = 189 GeV 



which shows a strong peak at 72.5 GeV owing to the 7Z-production subprocess. For the 
CC/NC process e^e~ —>■ uudd7, the electroweak diagrams dominate below the WW 
threshold, whereas the gluon-exchange diagrams dominate at the 7Z peak and above. 
The interference between purely electroweak and gluon-exchange diagrams is generally 
small. 

5 Summary and outlook 

The class of processes e+e~ —>■ 4 fermions + 7 has been discussed in detail. After 
classifying the different final states according to their production mechanism, the sets of 
all Feynman graphs are reduced to two generic subsets that are related to the two basic 
graphs of the non-radiative processes e"'"e~ 4 fermions. In this way, all helicity matrix 
elements are expressed in terms of two generic functions. Using the Weyl-van der Waerden 
spinor formalism, we have given compact expressions for these functions. We wrote two 
independent Monte Carlo programs, both using the multi-channel integration technique 
and an adaptive weight optimization procedure to reduce the Monte Carlo error. The 
results obtained with the two Monte Carlo programs agree well within the integration 
error. 

The detailed discussion of numerical results comprises a survey of integrated cross 
sections and photon-energy spectra for all different final states. Moreover, we have nu- 
merically compared different ways to introduce finite decay widths of the massive gauge 
bosons. Similar to the known results for e"'"e~ 4 fermions, we find that the application 
of running gauge-boson widths leads to totally wrong results for the radiative processes. 
Using constant widths consistently, leads to meaningful predictions. For the considered 
observables, the results for constant widths practically coincide with those obtained in a 
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e+(p+) 




W{q2) 



Z/e(/C2) 



e 



Figure 6: An example for a multi-peripheral diagram contributing to the process e"^e —>■ 

complex-mass scheme that fully respects gauge invariance. In the latter scheme gauge- 
boson masses are treated as complex parameters everywhere, in particular, leading to 
complex couplings. Similar to the situation for the well-known non-radiative processes, 
we find that for precise predictions the diagrams corresponding to triple-gauge-boson- 
production subprocesses are not sufficient and the inclusion of the other graphs is manda- 
tory. Finally, we have investigated the relevance of gluon-exchange contributions. In 
general, both the purely electroweak contributions and the gluon-exchange contributions 
are relevant and either one can be dominant depending on the process and the considered 
observable. The interference of both contributions is at most at the per-cent level. 

In this paper, we have assumed that the radiated photon appears as a detectable 
particle in the final state, i.e. soft and coUinear photons are excluded by cuts. The 
inclusion of soft and coUinear photons is, however, necessary if e"'"e~ 4fermions -|- 7 
is considered as a correction to four-fermion production. The analytical results of this 
paper and the constructed Monte Carlo programs will be used as a building block in the 
evaluation of four-fermion production in e"'"e~ collisions including 0{a) corrections. 

Appendix 

An example for phase-space generation 

To illustrate the phase-space generation in one channel, we choose the multi-peripheral 
diagram, shown in Figure ^ We determine the momentum configuration ki from the 
phase-space variables Sj, tj, (pi, and 6i with a suitable choice of mappings and calculate 
the resulting probability density g. 

First we study the propagator structure of the Feynman diagram of Figure |^ and 
choose a set $ of 11 phase-space variables that is suitable for importance sampling. Our 
choice of $ consists of three time-like invariants 



■S345 = ^345 = (^3 + ^4 + k^Y, 
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two space-like invariants 



^1 = ql = (p- - 



(A.2) 



and six azimuthal and polar angles 0i,2,3,4 and ^3^4, respectively. As worked out in detail 
in Ref. [ 3^, the 11-dimensional phase-space integral in (3J.) can be written as follows: 



.1=1 



5(4) 



V++V- 



1=1 J 



dS35 dS345 dSi345 R-?.{s) Rois 



V 



2 1 ■'1345, 



R2{sMb) R2{S35), 



(A.3) 



where -y/i denotes the CM energy with s = (p_ + The phase space of the 2 — > 5 

scattering process is composed of the phase spaces of two 2^2 scattering processes with 
two subsequent 1—^2 decays described by the following functions: 

R2{s): First, we consider the 2 — > 2 scattering process p+ +p_ — > A;i345 + ^2 in the system 
p+ + P- = where the direction of p+ is chosen to be the positive z axis. 

R2 (-51345): Then, we boost and rotate to the system A;/345 = and describe the 2-^2 
scattering process p+ + — > + k[, where the direction of p\ is chosen to be 
the positive z axis. 

R2( ■S345): Next, we define the 1 — >■ 2 decay of the Z boson /cg^g — > k'^^ + k'l in its rest frame 



h" — n 

^^345 ~ 



-R2(s35): Finally, we treat the subsequent 1-^2 decay of the virtual muon fcgg 0.3 
in its rest frame /C35 = 0. 

For massless external fermions, R5{s) explicitly reads^ 



h'" -I- h'" 



R^is) 



ds-. 



ds-. 



35 / ^1-5345 

35, cut "'S35 •^S34J 



ds 



1345 



As 



dti / d0i 
Jo 



X 



X 



4AV^(Sl345,tl,0) Jt2.n,in ^0 

A^/^(3345,^35,0) , , • 

/ dcos6'3 / d03 

0S345 J-l JO 



1 



1 j-2-K 

dcos^4 / d04 



-1 







with the kinematical function 

A(x, ?/, z) = + + z^ — 2xy — 2yz — 2zx. 
The limits on ti 2 are the physical boundaries of the 2 — > 2 scattering processes 

^l.min ~ -51345 ^l,max ~ 0, 



(A.4) 



(A.5) 



^2,min — (^1 " -51345) (■§1345 ~ -5345)/-5l345; t 



2, max 



0. 



■^As usual we set / = outside the imposed cuts. However, to improve the efRciency of the phase-space 
generation, cuts can be included already in the physical boundaries of phase space. In our example we 
have introduced a lower cut on S35. 
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In our example, we only have to cope with massive space-like propagators, so that we gen- 
erate the invariants S35 and S345 according to ( p.5|) and (|3.6|) , but do not apply importance 
sampling in the variables S1345, ti^2, 0i, 2,3,4, and cos 6*3 4, which are generated according to 
The resulting density g of (|3.9|) for this specific channel reads for z/ 7^ 1, 



■^^^^(■5345; -535, 0) 



g 42 82sS345 Al/2(si345,ti,0) 

S'-'' - 4'cut [(^345 - M|)2 + M|r|] 

X (S - S345) |il,mm| |i2,mm| 4 (27r)^ (A.6) 

with 1/1,2 = arctan[(si,2 - M|)/(Mzrz)], ■5i = S35, and S2 = s. Finally, the event charac- 
terized by the four- momentum configuration ki,i = 1, . . . , 5, and the beam momenta are 
defined in terms of the generated invariants and angles as follows: 

p^; = ^(1,0,0,1), 

= ^(1,0,0,-1), 

^2 = - — ^^^^^^^^-^(1, — COS 01 sin ^1, — sin 01 sin 6*1, — COS ^i) (A. 7) 

2^/ s 

with cos^i = (ti + 2p° A;0)/(2p° A;°), 

J >fi -^^^^(■51345 5-5345,0) A ■ n ■ a ■ n n \ / \ Q\ 

ki = —= (1, COS02 sm 6'2, sm02 sm6'2, cos6'2) (A-o) 

Sl345 



with COS 02 = it2 + 2p'^k[^)/{2p'^k[^) and p'^ = (^1345 - h)/{2^ 



1345 J 



k't = - — -COS03 sin 6^3, - sin 03 sin 6^3, - cos 6^3), (A. 9) 

^V"^345 

and 

^w/i _ v^^(^x, COS04 sin 6^4, sin04 sin6'4, COS04), 

fcr= (A.IO) 

The four-momenta fci, fcs, A;4, and k^ in the CM system are obtained from the four- 
momenta fci, fcg", k'l, and k'l' by applying the appropriate Lorentz transformations, as for 
instance described in Ref. [ ^ . 
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